clc,clear,close all
load("map_tokyo.mat")
%% 已知参数
Humidity = 0.59;
hight = y_altitude;
theta = 0;
V_wind = 3.*ones(1,length(GR));
Temperature = 288.75;

% figure()
% plot(GR)
%% TimeSpecialist man
CP = 411.1;
W = 20720;
s=0:250:750;
for n=1:length(s)
    tensity = s(n);
[P_plan,P_real,V,delta_T,fatigue,t(1,n)] = get_P(GR,CP,W,tensity,hight,V_wind,theta,Humidity,Temperature,turn_sign);
%需要画图可以去掉相应注释
% title('GR-x')
% figure()
% plot(V)
% title('V-x')
% figure()
% plot(P_real)
% title('P-x')
end

figure()
plot(s,t,'o-')
xlabel('sigma')
ylabel('t')
title('TimeSpecialist的时间t随sigma的变化')
%% Rouleur man
CP = 367.2;
W = 24640;
s=0:250:750;
for n=1:length(s)
    tensity = s(n);
[P_plan,P_real,V,delta_T,fatigue,t(1,n)] = get_P(GR,CP,W,tensity,hight,V_wind,theta,Humidity,Temperature,turn_sign);
%需要画图可以去掉相应注释
% title('GR-x')
% figure()
% plot(V)
% title('V-x')
% figure()
% plot(P_real)
% title('P-x')
end

figure()
plot(s,t,'o-')
xlabel('sigma')
ylabel('t')
title('Rouleurt的时间t随sigma的变化')
